arXiv:1505.07576v2 [math.AP] 20 Oct 2015 


STABILITY OF AN EULER-BERNOULLI BEAM WITH A NONLINEAR 
DYNAMIC FEEDBACK SYSTEM 
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Abstract. This paper is concerned with the stability analysis of a lossless Euler-Bernoulli 
beam that carries a tip payload which is coupled to a nonlinear dynamic feedback system. 
This setup comprises nonlinear dynamic boundary controllers satisfying the nonlinear KYP 
lemma as well as the interaction with a nonlinear passive environment. Global-in-time well- 
posedness and asymptotic stability is rigorously proven for the resulting closed-loop PDE- 
ODE system. The analysis is based on semigroup theory for the corresponding first order 
evolution problem. For the large-time analysis, precompactness of the trajectories is shown 
by deriving uniform-in-time bounds on the solution and its time derivatives. 


1. Introduction 

Let us consider a linear homogeneous Euler-Bernoulli beam, clamped at one end and with 
tip mass at the other free end. The state of the beam at time t is described by its transverse 
deflection u{t, x) from the zero-state, where x G [0, L] is the longitudinal coordinate of the 
beam, see Figure The well known PDE for the motion of the beam reads as 

puu{t,x) +Kv}^{t,x) =Q, (1) 

with the mass per unit length p and the flexural rigidity A. The boundary conditions for the 
clamped end at a; = 0 are given by 

u(t,0) = M'(t,0) = 0, (2) 

and for the free end at x = L, we have 

L) + Pm"{t, L) = —Te (3a) 

Muttit, L) - Au''\t, L) = -/e, (3b) 

where J and M denote the mass moment of inertia and the mass of the tip mass, respectively, 
and —Te and —fe describe the external torque and force acting on the tip mass. Here and in 
the following, the notation Ut is used for the derivative with respect to the time variable t, and 
u' for the cc-derivative. 

In literature, there exists a number of contributions dealing with the design of boundary 
controllers to stabilize this type of system. To mention but a few, in [1] the asymptotic stability 
was shown using semigroup formulation and applying the La Salle Invariance principle. To 
obtain stronger, exponential stability, frequency domain criteria [2], Riesz basis property [3], 
g] or energy multiplier methods [5], g] were employed. In contrast to these works, which are 
mainly based on linear static and dynamic boundary controllers, this paper is concerned with 
the interaction of the Euler-Bernoulli beam 0-(i with a finite-dimensional nonlinear dynamic 
system. In particular, it is assumed that this system generates a reaction torque Te = Tep +Tep 
and a reaction force fe = /e,i + /e, 2 ) respectively. The reaction torque and force is composed 
of the response of a nonlinear spring-damper system 

Te,i = di{u^{t,L)) + ki{u'{t,L)) 
fe.i = d 2 (ut(t, L)) + k 2 {u{t, L)) 


1 


(4a) 

(4b) 
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and the response of a finite-dimensional nonlinear system with state Zj G , j = 1, 2, 

{zi)t = ai{zi)+bi{zi)u[{t,L) (5a) 

Te,2=Ci(zi) (5b) 

and 


{Z 2 )t = 02 ( 22 ) + b 2 {z 2 )ut{t, L) (6a) 

/e,2= 02 ( 22 ), (6b) 

which constitutes a strictly passive map from the time derivative of the tip angle u^{t,L) to 
the reaction torque Te _2 and from the velocity of the tip position Ut(t, L) to the reaction force 
/e, 2 , respectively. The functions aj, bj, cj, dj, and kj, j = 1,2 as well as their mathematical 
properties will be specified in detail in the next section. 

The motivation for the setup- ([^ is as follows: In literature, when designing a boundary 
controller for the system Q - ([^ it is usually assumed that the external torque and force 

/e directly serve as control inputs. In this case, it is well known that the system 0 - 0 

can be stabilized (even exponentially) by a simple (strictly) positive linear static feedback, see, 
e.g., H, 0. However, in real practical applications the external torque Tg and force fe must be 
generated by some (electromagnetic, hydraulic or pneumatic) actuators whose dynamics cannot 
be neglected in general. In contrast to the usual approach in literature, it is therefore assumed 
in this work that these actuators are not ideally controlled, meaning that they are not serving 
as ideal torque and force sources, respectively, but that they are controlled in such a way that 
the subordinate closed-loop systems of the actuators comprising the actuator dynamics and 
a corresponding feedback controller constitute finite-dimensional passive dynamical systems 
according to ([^ and 0. In summary, the system 0-0 may be interpreted as a feedback 
interconnected system with the lossless Euler Bernoulli beam 0 - 0 in the forward path 
and the passive spring-damper system 0 as well as the strictly passive system 0, 0 in 
the feedback path, see Figure]^ It is well known that the feedback interconnection of passive 
systems preserves the passivity, see, e.g., 0. This fact is often exploited in the controller design, 
see, e.g., [TO] . m. for the finite-dimensional case. However, in the infinite-dimensional case 
the analysis is typically confined to linear systems, see, e.g., m, m, [la, or very recently |I4j . 
Thus, with this work we want to take a first step towards an extension of the state of the art 
to the nonlinear case by still considering a linear PDE but allowing for a nonlinear ODE at the 
boundary. 



u'{t, L) 


Figure I. Euler-Bernoulli beam with tip mass. 

The goal of this paper is to prove the global-in-time wellposedness and, most of all, the 
asymptotic stability of the feedback interconnected system 0-0 according to Figure]^ For 
both aspects, we have to deviate from the strategy employed in the analogous linear model (in¬ 
troduced and analyzed in [HIE]) : In the linear case, the generator of the evolution semigroup 
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environment + controller 


Figure 2. Interconnection of the Euler-Bernoulli beam system to a passive 
spring-damper system and a strictly passive system in the feedback path. 


is dissipative, which readily yields large-time solutions. The nonlinear semigroup for 0-® 
is not dissipative (in the sense of m)- Hence, standard semigroup theory will first only yield 
local-in-time solutions, and the construction of an appropriate Lyapunov functional for Q - 
Q then shows their global existence. 

Asymptotic stability of the linear counterpart model is based on the precompactness of the 
trajectories, which can be obtained from the compactness of the resolvent for the generator. For 
nonlinear evolution equations there exist different criteria for the precompactness of trajectories: 
They all split the generator of the nonlinear semigroup into the sum of a linear part C and a 
nonlinear part JV. In m A has to be maximal dissipative, and Af has to be integrable along 
the solutions. See [TH] for an application of this result, and also Section in this article 
below. Another approach is due to m, where only local integrability of Af along trajectories 
is needed, however, the semigroup generated by C needs to be compact. Finally, in it is 
shown that the trajectories are precompact if the nonlinearity A/” is compact, and the semigroup 
generated by C is exponentially stable. Furthermore, we refer to [21] for further results regarding 
precompactness of trajectories. Unfortunately, none of the above results apply to the problem 
discussed in this paper, except for the special case kj = 0 discussed in Section The reason 
for this is that the semigroup generated by the linear part is neither exponentially stable nor 
compact, and the nonlinearity A/” is generally not integrable along solutions. Hence, for Q - Q, 
we shall follow a different strategy, which was devised for a simpler system in [22] (it consists 
of a beam with a nonlinear spring and damper at the free end). For the precompactness of the 
trajectories of 0-0. we shall here prove uniform C^-bounds (w.r.t. time) on the solution, 
combined with compact Sobolev embeddings. 

Note that the beam in 0-0 (and in its linear counterpart) is undamped. Damping of 
the complete feedback system is only introduced via the damper of Q and the strictly passive 
systems 0, 0. This motivates that the linear model from m is asymptotically stable, but 
not exponentially stable. Hence, exponential stability also cannot be expected for our nonlinear 
system 0-0. Of course, exponential stability could be enforced by including damping terms 
into 0 (either a viscous damping of the form +aut or a Kelvin-Voigt damping of the form 
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+aU(^). While viscous damping would lead to a simple extension of the subsequent analy¬ 
sis, the higher order derivatives in the Kelvin-Voigt damping would require a rather different 
mathematical setup. Hence, we shall not elaborate on such dampings here. 

This paper is organized as follows: In Section the technical assumptions made for the 
coefficients and functions of the system 0 - 0 are specified, and in Section the problem 
is formulated as a first order evolution equation. Using semigroup theory we prove in Section 
1^ that it has a unique global-in-time solution. In Section the possible w-limit set of this 
evolution is derived and analyzed. For proving the asymptotic stability of ( [T|) - Q, we have to 
distinguish between two cases. For linear functions kj, it is shown in Section!^ that asymptotic 
stability can be achieved for all mild solutions. For nonlinear kj, it is much more involved to 
prove precompactness of the trajectories. In this case, asymptotic stability of classical solutions 
is shown in Section [T] 


2. Preliminaries 


In the following sections, we will give a rigorous mathematical analysis of the feedback inter¬ 
connected system 0-0 according to Figure]^ For this, the assumptions on the parameters 
and functions appearing in 0-0 have to be specihed. First of all, let us assume that the 
mass per unit length p, the flexural rigidity A, the mass moment of inertia J, and the mass 
M of the tip mass are constant and positive. For the spring-damper system Q we make the 
following assumptions for j = 1,2: 

(Al) There holds dj G (7^(11; K) (i.e. the space of two times continuously differentiable 
functions, see [23]), anc0 


d,(0) = 0. 

(7a) 

dj{s)s >0, Vs G M, 

(7b) 

d'j{0) > 0. 

(7c) 

Note that this implies dj{s)s > 0 for all s ^ 0. 

(A2) We require kj G (^^(K;]^), with kj{0) > 0 and 


Vk.{s) := f kj{a)d(7>0, VseK\{0}. 

(8) 


Based on the assumptions (Al) - (A2), it can be easily shown that the spring-damper system Q 
is strictly passive from the inputs u((t, L) and Ut{t, L) to the outputs Tep and /e,i, respectively, 
with the positive definite storage functions Vk-, j = 1,2, according to (j^. 

As a further consequence, we find uniquely determined constants Dj,Kj > 0 and functions 
Sj,Kj G C^(]R;K) with (5j{s) = 0{s^) and Kj{s) = O(s^) for s —0 such that 


dj(s) = DjS + Sj(s), Vs G K, (9) 

kj(s) = KjS + Kj(s), VsGR. (10) 


Hence, Djs is the linearization of dj{s), and KjS is the linearization of kj(s) around s = 0. 

(A3) Furthermore, we assume that there there exist (storage) functions Vj G 
for j = 1,2, such that for all Zj G : 

V,(0)=0, v,(z,)>0, (z,^0), 

VV,- (zj)- aj(zj ) < 0, (zj ^0), 

• b]{zj) = Cjizj). 


^Here, and k'^ denote the derivative with respect to the variable s. 


( 11 ) 

( 12 ) 

( 13 ) 
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According to the Kalman-Yakubovich-Popov (KYP) lemma for nonlinear systems with affine 
input, see Lemma 4.4 in [53], this implies the strict passivity of the systems ([ 5 ]) and ([^. 

For the mathematical analysis we furthermore require for j = 1,2: 

(A4) Assume there holds 

(14) 


lim Vj{zj)=oo, 

\zj\^oo 

Pj :=Hess(yj)(0) > 0. 


(A5) Uj, bj ^ C 


'2/mni. TUn, 


) and 


aj(0) = 0, 

det Aj 7 ^ 0, 


(15) 

(16a) 

(16b) 


where Aj = Ja (0) is the Jacobian of aj at Zj = 0. 

2 m", ;K)^ 


(A6) Cj e C' 


and 


Cj{0) = 0. 


The assumptions (A3) - (A6) have the following implications, for j = 1, 2: 

• There exists a unique regular matrix Aj G and a function aj G (IR"'^J ) 


such that for all Zj G 


^j(zj) = AjZj + aj(zj) 


J '^3 ^ '-^3 \^3 > 1 


|aj(^i)l = asZj-^0, 


(17a) 

(17b) 


hence AjZj is the linearization of aj{zj) around the origin. By using the first order 
Taylor expansion of VV} around the origin we conclude from (12) and (111 that 


^3 (^ 3 ^ 3)^3 — '^^3 G ® , 


and from (151 and (16) we find 


|VV;,-(zj) • aj{zj)\ > Clzjl"^ as Zj ->■ 0, 


(18) 


(19) 


for some positive constant C. 

• There exists a unique vector Bj G and a function Pj G (M”^; ) such that for 


all Zj G 


bjizj) = Bj + Pjizj), 

l/3j(2:j)l = Cld^il) asz^-^0. 

• There exists a unique vector Cj G and a function jj G C 


2/TIjn,- 


(20a) 

(20b) 

such that for all 


~ ^3 ’ ^3 ”1” 7i(^t), 

l7i(^j)l = C>(kir) aszj^O. 


Note that (13) implies 


(21a) 
(21b) 

(21c) 

To illustrate that the above assumptions may be satisfied we just mention the linear model 
from naini with aj = pj = 7 j = 5j = Kj = 0. There (and in many nonlinear perturbations of 
it) assumptions (Al) - (A6) hold. 


P 3 B 3 = C,. 


^Note that condition can be relaxed, see Remark 


4.10 
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Remark 2.1. For this paper it would even be possible to only make the weaker assumptions 
aj,bj e Cj G anddj,kj G Wjl’^(R;R), for j = 1,2. Here, 

is the space of all C^-functions whose first order derivatives are locally Lipschitz continuous 
functions (cf. |23jl. In particular, the local behavior of the functions aj, /3j, jj, 6j and Kj 
around the origin also stays the same, which can be seen by using the integral form of the 
remainder in the respective Taylor expansions. 

For the rest of this paper the conditions (Al) - (A6) on the system 0-® are tacitly always 
assumed to hold. 


3. Formulation as an Evolution Equation 


System 0-0 is reformulated as an evolution equation in the Hilbert space 
n = {y= [u,v,zi,Z2,^,'il}r ■ u G Hq{0,L),v G L‘^{0,L),Zj G G 


where for n > 2 

HSiO, L) := {/ G iJ"(0, L) : /(O) = /'(0) = 0}. 

Note that we impose the function value and its first derivative only at the left boundary, 
i.e. X = 11. Hence, Hq differs from the standard Sobolev spaces iJ" and Hg. We refer to [53] for 
the Lebesgue space A^(/) and the Sobolev spaces H^{I) on some interval I. The inner product 
is defined by 


{y, y)'H=-^ / dx + p 


1 ~ 1 ~ 
vvdx + + —i/u/> 


+ Kiu {L)u {L) + K 2 u{L)u{L) + zjPiZi + zjP 2 Z 2 , 


( 22 ) 


where the positive definite matrices Pj, for j = 1,2, are due to (15). 


operator 


For the following, the 


u 


V 

V 


P 

Zi 


ai{zi) + ^bi{zi)^ 

02(^2) + ]g& 2 (^ 2 )V’ 

Z2 

1— y 

e 


-Ku"{L) - [ci(zi) + di(|) + ki(u'{L))] 

_ Ku'"{L) - [C2(Z2) + d2(|) + k2{u{L))] _ 




is introduced on the domain 


D{A) = {yGn:uG H^{0,L),vG Ro(0,L), (23) 

^ = Jv\L),'ip = Mv{L)}. 


Based on the formulation of the coefficient functions, the operator A is decomposed into a linear 
and a nonlinear part: 

Linear part: The linear part is denoted by C, which is the linearization of A around the 
origin: 


u 


V 

V 


p , 

Zi 

1 _^ 

Ai^i + 

A2Z2 + 

Z 2 

1 7 

e 


-ku''{L) - [Cizi + + KW{L)] 



_Ku"'{L) - [C2Z2 + ijD2i} + K2 u{L)]_ 


and the domain is D{£) = D{A). 
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Nonlinear part: The nonlinear part JV is defined as the following continuous operator on 
all of H: 


u 


0 

V 


0 

Zi 


<Ti(^i) - 1 - 7/3i(^i)^ 

0 : 2 ( 22 ) + m/32(22)V' 

Z2 



-7i(zi)-(5i(5)-«:i(u'(£)) 



- 72 ( 22 ) - < 52 (m) - 


On D{A) there holds A = C + Af. 


Theorem 3.1. The linear operator C with domain D{C) generates a Co-semigroup of contrac¬ 
tions in TL, denoted by {e*^)t>o- 

Proof. This result has been shown for the same operator in Section 4.2 of m- For convenience 
of the reader we briefly sketch the main steps of the proof. A brief calculation yields for 


using (18): 


{Cy,y)'H =zj {PiAi)zi + zj (£ 2 ^ 2)22 

-D,\vfL)\^-D2\v{L)\^ <0. 


Hence the operator C is dissipative in P with respect to the inner product (22 1 . Furthermore, 
the inverse exists and is bounded (even compact). Now the statement immediately follows 
from the Lumer-Phillips theorem. □ 

Remark 3.2. Since C is the infinitesimal generator of a Co-semigroup of contractions, C is 
dissipative and ran(A — L) = P for all A > 0. In particular ran(/ — C) = P. So £ is hyper- 
dissipative according to Definition 2.1 in m- And Theorem 2.2 in m shows that £ is maximal 
dissipative, i.e. £ is not contained in a strictly larger dissipative operator (in the sense of graphs). 
This property is needed for the proof of Theorem |6.3| 


4. Existence of Solutions 

We are interested in solutions of the following initial value problem in P: 

yt{t) = Ay{t) = Cy{t)+Ny{t), 

2/(0) = yo&P- 


(24a) 

(24b) 


Any (mild) solution y{t) G C([0,r];'H), for T > 0, is known to satisfy Duhamel’s formula: 


y{t) = e‘^ 2/0 + f eO ^'>^J\fy(s) ds, 0 <t <T. 
Jo 


(25) 


Proposition 4.1. For every yo € P, there exists some maximal 0 < rmax(2/o) ^ 00 such that 
(24) has a unique mild solution y(t) on [0, Tniax( 2 /o))- If Ho £ D{A), the corresponding mild 
solution y{t) is a classical solution. //Tinax( 2 /o) < 00 , then Il 2 /(^)llw = oo- 

Proof. By assumption, the functions aj , jdj , 7 ^, 6j and Kj are continuously differentiable and 
locally lipschitz continuous, so the nonlinear map Af : P ^ P has the same properties. Fur¬ 
thermore, £ is the generator of a Co-semigroup, see Theorem |3.1| Now we may apply Theo¬ 
rem 6.1.4 in [25] to the autonomous problem (24), which yields the existence of a unique mild 
solution on the maximal time interval [ 0 , Tinax( 2 /o))- If Tmaxiuo) < 00 , then a blowup of yft) 
occurs. Moreover, Theorem 6.1.5 in implies that for yo G D{A), any mild solution is a 
classical solution. □ 










STABILITY OF AN EBB WITH A NONLINEAR DYNAMIC FEEDBACK SYSTEM 


Next we introduce the functional H : H ^ M., given by 

1 


(A|m'T + dx + + 


2M 


b(L) 

fci(s) ds - 


u(L) 

k2{s) ds■ 


Vlizi)+V2{z2). 


Vk^iu'iL)) Vk^{u{L)) 

Note that the first integral term in H(y) corresponds to the strain energy and kinetic energy of 
the Euler-Bernoulli beam, the next two summands are the translational and rotational part of 
the kinetic energy of the tip mass, 14^, for j = 1, 2, is the potential energy stored in the nonlinear 
spring elements, see Q and and Vj, for j = 1,2, are the non-negative storage functions 
of the strictly passive systems ([^ and ([^, respectively. Obviously H{y) > 0 for all y G V.. 
Note that H{y) is exactly the sum of the storage functions of the lossless Euler-Bernoulli beam 
(§-(§, the nonlinear spring-damper systems and the strictly passive nonlinear dynamic 
feedback systems ([^ and Q, cf. Figure}^ In the following, it will be shown that H qualifies 
as a Lyapunov function for the system ( |24[ ). 

Lemma 4.2. The function H is continuous in TL. 

Proof. The continuity of the terms in H except for the fcj-terms is immediate. Due to the 
continuous embedding H^{0,L) ^ L]) the continuity of the remaining fc^-terms follows 

as well. □ 


Lemma 4.3. Due to assumption (14) we have for any sequence {?/ra}neN C P: 


sup H{yn) < oo 

nGN 




sup WynW-H < OO. 

new 


Proof. It suffices to notice that {l^(zy„)}„gN is unbounded iff is unbounded. □ 


We now define the generalized time derivative along the mild solution y(t) of (24), i.e. for 
any initial value yo G H: 

t \,0 t 

which may take the value —oo. 

Lemma 4.4. For y^ G D{A) we have H{yo) = ^P[{y(t))\t=o < 0, i.e. FI is non-increasing 
along classical solutions. Here, denotes the right derivative. 


Proof. For yo G D{A) the corresponding solution y(t) of (24) is classical, and therefore has a 


continuous right derivative on [0, ri„ax(j/o))- So we can directly compute 

^(yo) = 

= ai(zi) • VVi(2;i) -I- 02(2:2) ■ VV2(2:2) - di{v\L))v’{L) - d2{v{L))v{L). 


Thereby we have used (13). The non-positivity of the generalized time derivative of the storage 
function H can be directly concluded from (12) and ([^. Clearly, this is also a consequence 
of the passivity property of the feedback interconnected system according to Figure This 
concludes the proof. □ 


Corollary 4.5. For y^ G D{A) the corresponding classical solution y{t) of (24) is global, i.e. it 
exists for all t G [0, 00 ). 


Proof. According to Lemma 4.4 H is non-increasing along y(t). Thus, according to Lemma 4.3 
no blowup occurs in y(t), and we have according to Proposition 4.1 that Tmax(2/o) = 00 . □ 
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Since Af is locally Lipschitz continuous and D{A) C H is dense, we can apply Proposi¬ 
tion 4.3.7 (ii) in [26] for the approximation of mild (non-classical) solutions: 

Proposition 4.6. Let i/q G H and {yo.njneN C D{A) be sueh that yo,ri j/o AL. Denote 
by yn{t) the global classical solution of (24) to the initial value j/o,n and by y(t) the mild 
solution corresponding to the initial value yo. Then yn{t) —t y{t) in C'([0,T];'H) for any T G 
(0,Tinax(yo))- 

Theorem 4.7. For any y^ G TL the corresponding solution y(t) of the initial value problem 


(24) is global in time. Furthermore, 1 1 —^ H{y(t)) is non-increasing on and y(t) is uniformly 
bounded in H on [0, oo). 

Proof. Consider yo G H and a sequence {yo,n} C D{A) with i/o,n —t Vo in AL. Due to the 
convergence yn{t) —t y{t) for all t € [0, T'max(2/o)) shown in Propositionand the continuity of 
F[, we get H{yn{t)) -G H{y{t)) for all 0 < t < rmax(yo)- Since H is non-increasing along every 
ynit), this implies also that 1 1 —)• F[{y{t)) is non-increasing on [0, Tniax(2/o))- Thus, according to 


Lemma 4.3 no blowup of y{t) can occur at t = T'max(2/o)- So, according to Proposition 4.1 the 
solution is global in time. Uniform boundedness of y{t) now follows from Lemma 4.3 □ 


Corollary 4.8. The function H is a Lyapunov function for the initial value problem (24). 


Remark 4.9. Since all mild solutions are global. Proposition |4.6| holds for any T G (0, oo). 

For every yo G FL and the corresponding mild solution y{t) we define S{t)yQ := y(t) for all 
t > 0. The family S = {S{t))t>o is a strongly continuous semigroup of nonlinear (bounded, 
continuous) operators in FL, cf. Theorem 9.3.2 in [26]. 


Remark 4.10. Since (14) is only needed to show that no blowup of the solution occurs, we 


may replace it by the weaker assumption 


lim Vj{zj) > H{yo), 

tj |—>-oo 


(14) 


depending on the initial condition yo for the problem (24). Thereby we argue as follows: 


According to Theorem 4.7 the function t i—> H{y{t)) is non-increasing (this is independent of 
(HI), which ensures that no blowup can occur in any component of y{t) except for Zj. If now 
Zi{t) or Z2{t) would blowup, we would get FF{y{t)) > H[yo) according to (fldj). So 

H{y{t)) could not be non-increasing, which is a contradiction. So (HI) is sufficient to show 
that no blowup occurs and that the solution is global in time. 


5. w-LiMiT Set 

In the following, S is the strongly continuous (nonlinear) semigroup generated by A on "H, 
defined at the end of the previous section. In this section, we investigate possible w-limit sets of 
S. However, non-emptiness of the w-limit sets will only be discussed in the subsequent sections. 
For yo G FL we define the trajectory 7 ( 1 / 0 ) by 

7 ( 2 / 0 ) := U S{t)yo. 

t>o 

Definition 5.1 (w-limit set). Given the semigroup S, the w-limit set for yo G FL is denoted by 
01 ( 2 / 0 ), and is the following set: 

^iVo) := {y GFL : 3{2„}„gN C M"'", lim tn = -koo A lim S{tn)yo = y}- 

n—FOQ n—Foo 

It is possible that 01 ( 2 / 0 ) = 0- 

According to Proposition 9.1.7 in [26] we have: 

Lemma 5.2. For y^ GFL the set 01 ( 2 / 0 ) is S-invariant, i.e. S{t)ui(yo) C o>(2/o) for all t > 0. 
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Let US consider now some fixed yi^ €%. According to the results of Section]^ the function 
t i-A H{S{t)yo) is non-increasing, and bounded from below by 0. Therefore, the following limit 
exists: 

^(j/o) := lim H{S{t)yo) > 0 . (26) 

t—¥00 

Lemma 5.3. Suppose oj{yo) ^ 0. Then there holds 

V2/Gw(i/o): H{y)=S]{yo). 

In particular, H{y) =0 for all y G uj{yo)- 

Proof. For every y G oj(yo) there exists a sequence {t„} C such that S{tn)yo -G y. Since 
H is continuous, cf. Section]^ this implies that H{y) = lim„_>oo-f^(5'(tn)2/o)- Due to (26) the 
right hand side equals i?(j/o)) and the result follows. □ 

We can use this lemma to identify the possible w-limit sets by investigating trajectories along 
which the Lyapunov function H is constant. 

Lemma 5.4. Let y GTL such that H{S{t)y) = f){y) for all t > 0, i.e. H is constant along 7 ( 1 /). 
Then 7 ( 1 /) C {y G TL : y = [u,v,0, 0,0, 0]^}. 


Proof. First, let y G D{A). We know from Lemma 4.4 and the corresponding proof that for all 
t > 0 


— H{S(t)y) = ai{zi) ■ Wi{zi) + 02 ( 2 : 2 ) • W 2 {z 2 ) 
- di{v'{L))v'{L) - d 2 {v{L))v{L), 


(27) 


where we omitted the dependence on t on the right hand side, i.e. \u, v, 21 , 22 , Jv'{L), Mv{L)]^ = 
S{t)y. Now (27) is required to be zero, and according to (12) and ([^ this holds iS ^ = ijj = 

2l = 22 = 0. 

Now let y G "H \ D{A). Then there is a sequence C D{A) such that yn ^ y as 

n -G 00 . According to Remark 4.9 we have S{f)yn —^ S(f)y uniformly on [0,T] for any T > 0. 
Therefore, we have also for the components 


Zj,nit) -G Zjit), 

in C'([0,T];M”0, 

(28) 

MVnit,L) -)• ipit), 

inC([0,T];K), 

(29) 

Jv'^it,L) -G fit), 

in C'([0,T];K). 

(30) 


Together with (27) this implies 




riGF 


is a Cauchy sequence in (^([O, T]; K). Since H is locally Lipschitz continuous in TL, we also have 
that {H{S{t)yn)}nGN is a Cauchy sequence in (7([0,T];K), so altogether {H{S{t)yn)}nGn is a 
Cauchy sequence in C'^([0,T];K). So there exists a unique h{t) G C^([0,T];]R) such that 


H{Sit)yn)^h{t) inCi([0,T];] 


(31) 


On the other hand we know that lim„_).oo H{S{t)yn) = H{S{t)y) = S){y) for every t > 0, and 
hence h{t) = Sj{y). Together with ( [^ this implies -^H{S{t)yn) —>■ 0 uniformly on [0,T]. By 
using (27) for every y„ this now yields that in ( [^ - (30) we obtain the limits Zj{t) = fit) = 
tpit) = 0. So Sit)y has to be of the form S'(t)y = [^(t), u(f), 0, 0, 0,0]^. □ 

Before we prove that the cu-limit set consists only of the zero solution, we need the following 
technical lemma: 
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Lemma 5.5. Let S be the nonlinear semigroup generated by A. For every yo G LL and for all 
t > 0 there holds: 


S{s)yods e D{A), 


and 


Sit)yo - yo = C f S{s)yods+ f AfS{s)yods. 
Jo Jo 


(32) 


(33) 


For the proof we only need the fact that C generates a Co-semigroup, and that Af is dif¬ 
ferentiable and locally Lipschitz continuous. Hence, the above result still holds true for more 


general operators C and Af, which satisfy the mentioned properties. The proof of Lemma 5.5 


is analogous to the proof of Lemma 5.4 in [^, see also m for a general version of this lemma. 

Theorem 5.6. Let ib ^ Ll G H be an S-invariant set sueh that H\q is constant. Then H = {0}. 
In particular, for any y^ G H either w(?/o) =% or uj{yo) = {0}. 


Proof. Take a hxed yo G LI, and let y{t) be the corresponding mild solution of (24). Clearly, 
7 ( 2 / 0 ) C LI, and according to Lemma 5.4 y{t) is of the form y{t) = [u(f), u(t), 0,0,0,0]^. 

Step 1 (linear system for u{t), v{t)): First we note that, according to (|32|), there holds for all 
t >“(h 

0= / 'if>{s)ds = M j v{s, L) ds = M{u{t, L) — uo{L)), 

Jo Jo 

0= / ^(s)ds = j(" f u(s,x)ds') 

Jo ^Jo ' 

= J{u'{t,L)-u'^{L)). 


Thus u{t,L) and u'[t,L) are constant in time. According to (33) the (projected) mild solu¬ 
tion yp{t) = [^(t), i;(t)]^ satishes the following system (i.e. the hrst, second, fifth, and sixth 
component of (33)): 

pt 

u{ ■' 


i{t) — uq = / u(s,ai)ds, 

Jo 


Kt) -Vo = - 


A 


u{s, x) ds 


0 = A 

0 = -A 




J u{s, x) ds^ 

( / u(s,a;)dsj -|- / Ki{u' {s, L)) ds, 
^Jo ' Jo 

^ J u{s, x) ds^ 

^ J u{s, x) ds^ 


(34a) 

(34b) 

(34c) 

(34d) 


K 2 {u{s, L)) ds. 


Mild solutions satisfy u G C(U~'~; 11^(0, L)). Hence, we can interchange the integration and 
differentiation in the last term of (34c). Using the fact that u'{t,L) is constant, we have (for 
u'oiL) 7 ^ 0 ): 

rt 


Ki{u'{s, L)) ds = tKi(uQ(L)) 

ki{u'o{L)) 


u'q{L) 


^ J u{s, x) ds^ 


x—L 
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Next we define the constants (since Kj{0) = 0): 

Ki{u'oiL)) 


Ki := Ki + 


K 2 := K 2 + 


'^iL) ■ 

K2{uo{L)) 

Uo{L) 


if Uq{L) 7 ^ 0 , else Ki := Ki, 
if uo{L) 7 ^ 0, else K 2 ■= ^^ 2 - 


(35) 


With this we may rewrite (341 as 


i{t) — uo = / v{s) ds, 

Jo 


}{t) -Vq = - 


A 


l{s) ds 


0 = A 
0 = -A 


u(s, x) ds 


x—L 


Ki 


^ J u{s, x) ds^ 


X—L 


^ J u{s, x) ds^ 


X—L 


+ K 2 / u(s,a;)ds 


X—L 


(36a) 

(36b) 

(36c) 

(36d) 


making this system linear. Thus, the projected vector yp{t) = [^(t), i;(t)]^ is the unique mild 
solution of 


{yp)t — J^pypi 

yp{Q) = [■uo,uo]’^, 


(37a) 

(37b) 


with the operator 


r ■ 

*^p • 


u 


V 



A^JV 

V 


\ 


The equations (36c) and (36d) are incorporated into the domain D{Cp). For further details on 
the operator Lp in the space T-Lp see the Appendix. 

Step 2 (proof of u{t,L) = u'{t,L) = 0): We now investigate solutions of the projected prob¬ 
lem (37) with the additional property that u(t,L) and u'(t,L) are constant in time. Since the 
semigroup e*^p is unitary in Hp, we know in particular that ||?;(t)||i 2 < C = -^\\yp{0)\\'Hp for 
alH > 0 (cf. (69)). Applying the norm to (36b) this yields 

rt 


sup 

t>o 


i(s) ds 


L2(o,L) 


< 00 . 


(38) 


Next we apply the following Gagliardo-Nirenberg inequalities (cf. [H]); which guarantee the 
existence of a C > 0 such that there holds for all t > 0: 


u{s) ds 


ds 


L“(0.L) 


< C 


il 


L°°(0,L) 


< c 


* \ IV 

u{s) ds 


L2(o,L) 

3 

8 

L2(o,L) 


u{s) ds 


/o 


i(s) ds 


L2(0,L) 


L2(o,L) 


(39) 


The first factor on the right hand side in both inequalities is uniformly bounded (with re¬ 
spect to t) due to (38). For the second factor we observe that, according to Theorem |4.7[ 
t I—>■ ||u(t)||L 2 (Q 2 ,) is uniformly bounded, and therefore t >—t |j tt(s) ds|jj;, 2 (o grows at most 
linearly. Hence, (39) implies that t >->• || m(s, L) ds||ioo(g_^) grows at most like ti and 
ll/o w'(s.^)ds||L “(o.L) most like ts as t ^ 00 . But this contradicts the fact that u{t,L) 

and u'{t,L) are constant, unless Uo{L) = Uq{L) = 0. This shows that u{t,L) = u'{t,L) = 0 for 
all t > 0 . 
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Step 3 (Holmgren’s Theorem): By iterated t-integration we shall now construct C^-solutions 
of (37a|), for which we can apply the Holmgren Uniqueness Theorem [521 Section 3.5]. So we 


define yi{t) = ui(t)]^ := yp{s)ds + Due to Theorem 1.2.4 in [25] and 


Lemma A.l we have yi{t) S D{Cp) for all t > 0. So yi is a classical solution of ( |37a ) to the 
initial condition 2/i(0) = Furthermore, because of u{t,L) = u'{t,L) = 0, again 

ui{t, L),u'i(t, L) are constant in time. Completely analogous to the previous step we can show 
again that Ui{t, L) = L) = 0. 

Next we shall construct solutions of higher regularity. We iterate the previous step and define 
recursively ?/„(t) = [u„(t),:= ds + Up^luo,which solves (37a) classically 

with the initial condition j/„(0) = £“”[uo,uo]^. Again we have Un{t,L) = u'^it^L) = 0. 
Furthermore, by definition we have on the one hand Cpyn{t) = yn-i{t). And on the other hand 
jO-p[un,Vn]^ = [u„, — A/pu]7]^, SO we Can show inductively that ?/„ S C'(K+; L) x 

i7g"(0, L)). Now we consider the solution for n > 2. It satisfies the following partial 
differential equation with boundary conditions: 

/■ \ IV 

[Unjtt = --U„ , 

P 

[un{0,x), (■u„)t(0,a;)]^ = Cp'^[uo,vo]^, 


■u„(t, 0 ) = 0 ) = 0 , 

Un{t, L) = ... = L) = 0. 


(40a) 

(40b) 

(40c) 

(40d) 


By using (40a), Un G C'(K+; .Hq"''''^( 0, L)), and the fact that {u„)t = G C(M+; 77q"(0, L)), 


we obtain the following properties for the mixed fourth order space-time derivatives of Un- 

<''GC(M+;i72»-2(o,L)), 

(u„)"'GC(M+;i7o'”-3(0,L)), 

(w„)" = e 

(u„);,, = G C(M+;772"-5(0,L)), 

Mtttt = e C(M+;#2'^-®(0,L)). 


So for n > 4, all mixed derivatives of of order four lie in C'(K“''; Ho{0, L)) C C(K+ x [0, L]). 
Thus Un{t,x) is a C"*-solution of (40). 

Now we can apply the Holmgren Uniqueness Theorem [521 Section 3.5] on the strip K+ x 
(0, L). Due to (40d) all partial derivatives up to order 3 of U 4 vanish on the line M"*" x {L}. 
Therefore, Holmgren’s Uniqueness Theorem implies that U 4 = 0 has to hold everywhere in this 
strip. (See also the proof of Lemma 3 in [1] for a similar result - but without a detailed proof.) 
Therefore C~'^[uq,Vq\^ = 0 has to hold, and since C~^ is injective, this yields [uo,Uo]^ = 0. 
Since yp{t) = e*^i’[uo,uo]''^, we conclude that u{t) = v{t) = 0 for all t > 0, and hence D = {0}. 
For the final statement of the theorem, let oj{yo) 7 ^ 0. Then, by Lemma 5.2 uj{yo) is S- 


invariant, and by Lemma 5.3 H is constant on CLi(un). Hence, by the first statement of Theo- 
reml^ w(j/o) = { 0 }. □ 


As a consequence we obtain convergence to zero for trajectories with u}{yo) 0: 
Corollary 5.7. If uj{yo) ^ 0 for some yo G H, then 

lim ||S'(t)?/o||w = 0. 

t—>-oo 
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Proof. If uj{yo) ^ 0 then there exists a sequence {tnjriGN with —?> cx) such that 

lim„_>.oo S{tn)yo = 0- Due to the continuity of the Lyapunov function H this implies that 

lim H{S{tn)yo) = 0. 

n—>-oo 

But since t i-A H{S{t)yo) is non-increasing, this implies that even 

lim H{S{t)yo) = 0. 

t—>-oo 

Due to the continuity of H this implies that \\S{t)yo\\'H —t 0 as f —oo. □ 


6. Asymptotic Stability - Linear fc,- 


In the case where the kj are linear we are able to show precompactness for all trajectories, 
even for the mild, non-classical solutions. This will yield that the w-limit set is always 

non-empty, and hence the asymptotic stability of the nonlinear semigroup will follow. 


Lemma 6.1. Let y^ S %, and y(t) be the corresponding mild solution of (24). For j = 1,2 let 
Kj =0. Then Afy(t) € 

Proof. First, let us assume that yo & D{A), so y{t) is a classical solution. We know from 
Theorem 4.7 that Hijjit)) is non-increasing. By integrating (27) with respect to time we obtain 


H{yiT))-H{yo) = 


-di 




J J J 


— do 


M 


(41) 


-I- ai{zi) ■ VTi(zi) -I- 02 (^ 2 ) • 'VV 2 {z 2 ) dt =: Irivo), 


where all terms on the right hand side include elements of the vector y(t), thus depend on t. If 
we let T —7^ 00 , we know that P[{y(T)) —>■ be. the limit exists and the integral looivo) is 

hnite. 


Now we consider yo G 'hi, and y(t) is the corresponding mild solution of (24). Let {i/o,n}nGN C 
D{A) be a sequence with yo,n —t Vo- According to Proposition |4.6| and Remark [4^ the cor¬ 
responding classical solutions y„(t) converge to y{t) in C'([0,T];77) for all T > 0. Therefore 
Iriyo.n) Iriyo), cf. (|4^. Due to continuity of i7, also i7(y„(r))-i7(yo.n) ^ H{y{T))-H{yo) 
as n —)■ 00 . Thus, ( [iT] ) also holds for mild solutions for any T > 0. Since P[{y{T)) S^(yo) G 
[0, i7(yo)] as T —>• 00 , the integral Loo(yo) is finite. 


Now we know that for any (mild) solution y(t) the integral /oo(?/o) from (41) is hnite. Since 
all the terms in the integrand of ( [dlj ) are non-positive, we conclude together with (19) and 0 
that 

z,(t),m,m€L^(R+). (42) 

Under the assumptions we made in Section for the functions occurring in the nonlinear 
operator AC, the properties ([4^ immediately imply My{t) G L^(M+;'H). □ 


Remark 6.2. To obtain Afy(t) G L^(]R+;'H) in the above proof, we used in (41) that the 


nonlinear damping functions dj include a non-vanishing linear part (i.e. Dj > 0). The same 


assumption will also be needed in Step 3 of the proof of Lemma |7.2| below. However, in the 
nonlinear spring-damper system of |22j . a locally quadratic growth of the damper law was 
sufficient. From a practical point of view, this is not restrictive at all. 


We note that (41) does not give any control on u{t,L) and u'{t,L) (in the sense of ([4^). 
Hence, the linearity assumption Kj = 0 was crucial for the above proof. 


Theorem 6.3. Let Kj = 0 for j = 1, 2. For any yo G TL there holds limt_>.oo S{t)yo = 0, i.e. the 
semigroup S is asymptotically stable. 
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Proof. Our aim is to apply a version of Theorem 4 in m- It states that if £ is a linear, maximal 
dissipative operator with (A — C)~^ is compact for some A > 0, and / S L^(IR+;H), then every 
mild solution of the Cauchy problem y{t) — Cy{t) + f{t) has a precompact trajectory. 

According to Remark |3.2| the linear part £ of Al is a maximal dissipative operator on T-L. 
As seen in the proof of Theorem ^ exists and is compact. Since C generates a Co¬ 
semigroup of c ontractions, (A — C)~^ exists and is compact for all A > 0. Finally, according 
to Lemma 6.1 we know that Ny{t) G L^(K+;'H) for y(t) := S{t)yo. Due to these facts, we 
can apply Theorem 4 in m with f{t) := Ny{t). This shows that the w-limit set w(yo) is 
non-empty. Thus, due to Corollary 5.7 and Theorem |5.6[ we conclude w(yo) = {0} and that 
the entire solution y[t) converges to zero. □ 


7. Asymptotic Stability - Nonlinear kj 

According to Corollary |5.7[ any trajectory with a non-empty w-limit set already is asymp¬ 
totically stable. Thus, in order to complete the discussion we show in this section that (at 
least) any classical trajectory possesses a non-empty w-limit. We do this by proving that every 
classical trajectory is precompact. To this end we follow a strategy introduced in [52]. We 
begin with the following preparatory result (which would be obvious for linear semigroups): 


Lemma 7.1. Let y(t) be a (mild) solution of (24) and let yo G D{A^) := {y G D{A) : Ay G 
D{A)}. Then y G C'^{[0,oo)]'H) and yt(t) G D{A) for all t>0. 


Proof. If we already knew that y G (7^([0, oo); 77), it would follow that y := yt satisfies 


0 

0 

a[{zi)zi + ^[I3[{zi)zi£, + I3i{zi)i] 
a'2{z2)z2 + ■ii[l32{z2)z2lp + l32iz2)^] 
-ii{zi)zi - k;(u'(L))u'(L) 

-12{Z2)Z2 - - k'2{u{L)){l{L) 


(43) 


4.1 


However, at this point we only know that y{t) G C'^([0,oo);'H), see Proposition 
by (43) we define the following functions for this fixed y{t) = [u, v, zi,Z 2 , C, '4’V(t)' 


Motivated 


Gi{t, Z) a'i{zi)C,i + -j{l3'i{zi)C,ii + /3i(zi)S], 
G2{t,Z) := a'2{z2)C,2 + ^[P2{z2)C,2'lp + P2{z2)'^], 
G3(CZ) := -7'i(^i)Ci - - ^^\{u'{L))U\L), 

G^{t,Z) := -7^(z2)C2 - k'2{u{L))U{L), 


where Z := [t/, P, ^i, C 2 j 2, G TL. Since y{t) is a classical solution, it follows from the 
regularity assumptions of the coefficients that t Gj{t, Z) lies in for all j = 1,..., 4. As a 
consequence the operator M ■. [0, T] x 77 —>■ 77 defined by 

A7(7, Z) := [0,0, Gi(7, Z), G2(7, Z),G^(t, Z), G^it, Z)]^, 

is Lipschitz continuous for any fixed T > 0, and linear in Z G 77. Now the linear, non- 
autonomous initial value problem 

Zt= TZ ^M[t,Z), 

Z(0) = Zq&'H, 


(44a) 

(44b) 
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is considered. According to Theorem 6.1.2 in there exists a unique global mild solution 


Z{t) of (44) for every Zq S H. If Zq G D{L) this solution is classical, see Theorem 6.1.5 in [5S]. 

Our next aim is to show that for the classical solution y(t) fixed in the beginning, the 
(continuous) function yt{t) is indeed a mild solution of (44) for Zq = Ay^. Since y{t) satisfies 
the Duhamel formula ([2^ and is differentiable, we obtain after differentiating with respect to t 


(45) 


ytit) = A^Cyo + ^ e(‘-*)^AA 2 /(s) ds. 

According to the proof of Corollary 4.2.5 in the following statement holds true 

e^*-^^^Afy{s) ds = e*^A/'j/o + / -^Afy{s) ds. 

Jo 


(46) 


Inserting (46) in (45) yields that yt (t) fulfills the Duhamel formula for (44). As a consequence 


yt{t) is the unique mild solution of (44) to the initial condition Zq = Ayo- Moreover, we know 
that this mild solution Z{t) = yt{t) is a classical solution of (44) if Ayo G D{A), i.e. yo G D[A?). 
Hence yt G C'^(M+;H) and y G C'^(M+;H). □ 


Lemma 7.2. The trajectory 7 ( 1 / 0 ) precompact in TL for yo G D{A^). Moreover, there exists 
a constant C > 0 such that 

||j/*(t)||„<C, Vt>0, (47) 

where C depends continuously on ||2/oll« l|yt(0)||«. 


Proof. In order to prove precompactness of the trajectory, it suffices to show that 


sup ||Ay(t)||w < 00 , 

t>0 


due to the compact embeddings H'^{0,L) iJ^(0, L) L^(0, L). However, this is 

equivalent to showing that yt is uniformly bounded in "H, since yt = Ay. Again, this is equivalent 
to 

^(2/t) = ^^ «)^dx+^(u'„(L))V y(utt(L))^ 

rAA) fUt{L) 

+ / fci(s)ds+ / k 2 {s)ds+ Vi{{zi)t)+ V 2 {{z 2 )t) 

Jo Jo 

being uniformly bounded. Since y{t) is a classical solution, we have the following equalities 


Ut{L) = 


M’ 


n'tiL) = |. 


According to Theorem 4.7 those terms are always uniformly bounded. Moreover, due to regu- 

we see from (5a) and (6a) that {zj)t G L°°(M+) 


4.7 


larity of the functions aj , bj and Theorem 
for j = 1,2. Therefore, the boundedness of H{yt 
functional 


is equivalent to the boundedness of the 


«)^dx+^(u'„(L))V y(utt(L)) 


Hence, our aim is to derive a system of equations satisfied by yt(t), and then to show that H{yt) 
is uniformly bounded. 

Step 1 (Time derivative of the system): According to Lemma|7.ll y{t) G (^^([0, 00 ); "H). Dif¬ 
ferentiating IFur with respect to time hence shows that yt is the classical solution of the 
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following system 


puttt + = 0, 

0 ) = 0 , 

AUt{t,L) + JUttt(t,L) + (Te)t(t) = 0, 

—Au”{t, L) + L) + {fe)t{t) = 0, 


where 

Te ■■= ci(zi) + di{u^{L)) + ki{u'{L)), 

fe ■■= 02 ( 22 ) + d 2 {ut{L)) + k 2 (u{L)). 


Therefore, from (49) it follows 


{T^)t = Vci{zi)-{zi)t+d[{u[{L))uu{L)+k[{u' {L))u^{L), 
{fe)t = yc2{z2)-{z2)t+d2iut{L))uuiL) + k2iu{L))ut{L), 


and from (5a) and (6a), we obtain 


{Zl)tt — [Jai{zi) + Wt(I/)Jbi(2l)](zi)t + bl{Zl)u'^^{L), 
{Z2)tt = [JaAz2) + Ut{L).Jb^{z2)\{z2)t + b2{z2)utt{L), 


(48a) 

(48b) 

(48c) 

(48d) 


(49) 


(50) 


(51a) 

(51b) 


where Ja^, J bj de note the Jacobian matrices of the functions aj and 6^-, respectively. Note that 
: it follows that Zj{.),ut{. , L) = m((. ,L) = ^G L^(M+) (cf. (42) for a similar 

e 


from Lemma 


4.4 


conclusion). Therefore (5a) and (6a) imply {zj)t 
Step 2 (Time derivative of H{yt)): We obtain 




UtttUtt dx + A u'Lu'I dx 


“b '^^ttt(A)utt(T) + Muttt{L)utt(,L) 

= Utt{L){Muttt{L)-Au”'{L)) 

+ ^tt(A)(Aw"(I/) + Ju[f.^{L)) 

= -UttiL){{Z2)JyC2{z2) + k'2{u{L))ut{L)) 

- n;*(L)((2i)7Vci(2i) + k[{u'{L))u[{L)) 

- d'2{ut{L)){uu{L)f - d'M{L)){u'M)\ 


(52) 


where we have performed partial integration in x twice, and then used ( |48[) and (50). Integrating 
(52) on the time interval [0,t], for some arbitrary t G ffi"*", we get with (T^ 

+ hit) + hit), (53) 


where 

hit) := - J^UttiL)(^izi)J\'ciizi) + k[iu'iL))utiL)^ ds, 

hit) ■=- J Wtt(A) (( 2 : 2)7 Vc2 (22)+ fci('«(A))wt(A)) 


ds. 
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Step 3 (Boundedness of Ji and 12 )- Next, we show uniform boundedness for each component 
of I 2 by using partial integration in time: 

- J Utt{L)k 2 {u(L))ut{L) ds = {ut{t, L)f k 2 (ut{t, L)) + ^ (ut(0, L)f k^iutiO, L)) 

+ - / ut{Lfk”{u{L)) ds <C, yt> 0. 

2 Jo 


Further, it holds that 
rt 


Utt{L){z 2 )J'VC 2 (z 2 ) ds = Ut{t, L){z 2 )t{t) ' VC 2 (z 2 (t)) “ (0, L) (Z 2 )t (0) ' VC2(z2(0)) 


- f Ut{L)[{z2)J[ileSs(c2)iz2)]iz2)t + {z2)Jt'yC2iz2)]ds. 
Jo 

Since C 2 G and 22 (0 £ it follows that 

f \utiL){z2)i [iless{c2)iz2)]{z2)t\ds < C f I(z2)tpds, 

Jo Jo 

and (with @) 

/ UtiL){z2)JtVC2iz2)ds = f Ut{L)[Ja^{z2){z2)t +Ut{L)Jb^{z2)iz2)t]^VC2iz2)ds 
Jo Jo 

+ [ VC 2 iz 2 y b 2 {z 2 )uttiL)utiL) ds 

Jo 

= / Ut{L)[Ja2{z2)iz2)t+Ut{L)Jb2{z2)iz2)t]^VC2iz2)ds 

Jo 

+ ]^VC2{z2{t)y b2{z2{t))ut{t,Lf 

- ^VC2(z2(0))^&2(22(0))ut(0, L)^ 


)1 


-II (uALfiz,) 

■ [db2(z2)^Vc2(z2) +Hess(c2)(22)&2 

<c[ \utiL)\'^ + |( 2 ; 2 )tpds 

+ ^'^C 2 {z 2 {t))^b 2 {z 2 {t))utit,Lf 
- ^Vc2(z2(0))^&2(2:2(0))ut(0,L)^. 


ds 


For the estimate of the second integral we have used the uniform boundedness of (^ 2 )*, see 
the discussion before Step 1 of this proof. The uniform boundedness of Ii follows analogously. 
Hence, H{yt(t)) is uniformly bounded in time. Furthermore, it can be seen that all the positive 
constants C appearing in the above calculations depend continuously on the initial conditions. 
This concludes the proof. □ 


In order to extend this result to all classical solutions, we need the following density argument. 

Lemma 7.3. For any y G D{A) there is a sequence {i/n}neN in D{A^) such that lim„_>.oo Vn = y 
and lim„_>oo Ay-n = Ay. 
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Proof. Let an arbitrary y G D{A) be fixed. Notice that it suffices to show that there exists a 
sequence {yn}nen with = [un Vn zin Z 2 n in D{A'^) such that lim„_j.oo 2/n = 2 / in the 

space H*{0, L) x L) x x x M x K. The set D{A^) = {y G D{A) : Ay G D{A)} is 
equivalent to 


vGHI(0,L), (54) 

u G Hl{Q, L) A = m''(0) = 0, (55) 

^ = Jv'{L), (56) 

'll; = Mv{L), (57) 

^u^{L) = Au"(L) + [ci(zi) + di ^ + ki{u'{L))], (58) 

^u^^iL) = -Au"'iL)+[c2{z2) + d2 [^)+k2 {u{L))]. (59) 


Since (7^(0,L) := {/ € C°°[0,L] : /^^^(O) = 0,V/q € No} is dense in Hq{0,L) (see Theorem 3.17 
in m), there exists a sequence {v„}„gN C Clf°{0,L) such that liin„_>oo iin = u in i7^(0, L). 
Also, Vn satisfies (54), for all n G N. Defining := Jv’.^[L) and ifn ■= Mvn{L) ensures that 
yn satisfies (56) and (57). Moreover, the Sobolev embedding H^{0,L) ^ C^[Q,L\ implies that 
lim„_>oo f.n = f. and lim„_>oo lAn = V' as well. Next, let zin := zi and Z 2 n ■= Z 2 for all n G N. 

Finally, the sequence {u„}nGN C C°°[0,L] will be constructed such that it„ satisfies (55), 
(58), and (59) for all n G N, and lim„_>oo Wn = m in H'^{0,L). To this end we introduce an 
auxiliary sequence of polynomial functions 


.— ^2,n^ “t“ ^3,n^ 4“ ^6,n^ 4“ ^7,n^ 4“ ^8,n^ 

4 “ 4 " h^Q j^iX 4 “ hw jiX , 


for all n G N, where h 2 ,„,..., hu n G K are to be determined. It immediately follows that 

h„(0) = h'jo) = C(0) = hl{Q) = 0. (60) 

Let /i 2 ,ra = 2 °^ and /i 3 _„ = “ , which is equivalent to 

K{Q) = u"{Q), C(0)=ri"'(0). (61) 

Further conditions are imposed on /i„: 

hW(L) = wW(L), fcG{0,l,2,3}. 

This can equivalently be written in terms of coefficient^ 


Ti = + hr^nL + hg^nL^ 4" hg^nL^ 4- hiQ^nL^ 4- hii^nL^, (62a) 

72 = 4- Ihr^nL + Shg^nL^ 4- 9hg nL^ 4- 10/lio,nL^ 4" 11/lll^riL®, (62b) 

73 = 6-/i6,ra 4- 7-hr^nL 4" 8-hg^nL^ 4" 9-hg^riL^ 4" lO-hio^nL"^ 4- 11-hii^nL^ (62c) 

74 = 4- 7-hr^nL 4- 8-hg^nL^ 4- 9-hg^nL^ 4- lO-hio^nL"^ 4- (62d) 


with 


u{L) u"(0) u'"{0) 

~ L6 2L4 6L3 ’ 

u"(L) u"(0) u"'(0) 

~ ~lA L4 Z^’ 


u'(L) u"(0) u"'(0) 

^ “ L5 L4 2L3 
u"'(L) u"'(0) 

~L3 


^The coefficient k- (the Pochhammer symbol, see [30] ) for k, I G N, I < k is defined by k- := k ■ (k — 1) ■ 
I + 1 ). 


■{k- 
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We further require that satisfies: 
AM 


P 

AJ 


{L) — —Au'''{L) + 02(22) +k 2 {u{L)) 

kni^) — Au" {L) + Ci( 2 i) + di + ki(u\L)) 


= ■ T5, 


:= re, 


where (631 and (64) are equivalent to: 

AML2 ~ + '^^h'j^nL + S-hg^nL^ + 9-hQ^nL^ + 10-hiQ nL^ + II -/111 


re 


AJL 


— 6^^6,n + + S-hs^nA^ + 9-hg + 10-hio,nA'^ + II -/111 


(63) 

(64) 

(65a) 

(65b) 


Such hn exists and is unique, due to the fact that linear system (62) and (^65| has strictly 


positive determinant. Consequently, ( [60| , ( [6l| ), and (62) imply that m — /i„ G 1^(0, L), for all 
n G N. Since C^{0,L) is dense in Mg(0,L), there exists a sequence {M„}nGN C C^{0,L) such 
that \\un — {u — h„)||//4 < i, Vn G N. Now defining := + /i„, gives lim„_,,oo = it in 

H'^{0,L). Obviously ii„ satisfies (55) for all n G N. Also, due to (63) and (64), ii„ satisfies (58) 
and (59), as well. Hence, the statement follows. □ 

Theorem 7.4. For all yg G D{A) the trajectory 7 ( 1 / 0 ) *5 precompact in H. 

Proof. Let yg G D{A) be chosen arbitrarily, and let {i/no}neN C D{A^) be an approximating 
sequence as in Lemma [7.3| Then there holds: 

lim Ayno = Ayg. (66) 

n—^oo 

For an arbitrary T > 0, and by applying Proposition |4.6| it follows that the approximating 
solutions yn{t) converge to y{t) in C{\f),T]]'H). Since yn{t) G C'^([0,oo);'H) and solves (24) for 
all n G N, (66) yields 

lim (i/„)t(0) = Ayg in U. 


(67) 


Hence, (47) and (67) imply that there exists a constant C > 0 such that for all n G N: 

sup||(yn)t(t)|lw < Cdlyollw, Pl/olk), 
i>0 

where the constant C does not depend on n. From here it follows that {yn)t is bounded in 
L°°(K+;'H). Hence, the Banach-Alaoglu Theorem (see Theorem 1.3.15 in m) implies that 
there exists w G L°°(M+;'H) and a subsequence {yn^lfeGN such that 

{ynjt^w in L°°(M+;77). 

For arbitrary z G H and t > 0 there holds 


lim / ((j/„Jt(T), 2 )«dT = / (w(r), 2 )«dT, 
k->^°°Jg Jg 


which is equivalent to 

lim (y„fc (t) - yuk (0), z)u = { / w(t) dr. 
Since lim„_>oo Unij) = yij) (in TP) for all r G [0, 00 ), it follows that 

{ylt) - 1/(0), z)-H = (^J^ w{t) dr, 

Since z G H was arbitrary, we obtain 

yit) - 2/(0) = [ dr, Vt > 0. 

Jo 


n 


( 68 ) 
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Due to continuous differentiability of y, the time derivative of (68) can be taken, which yields 
yt = w. This implies yt £ i.e. ||2/t(-)ll -H is uniformly bounded, which proves the 

theorem. □ 


Corollary 7.5. For any yo £ D{A) there holds limt_>.oo = S{t)yo = 0. 


8. Conclusions 

In this paper, we provide a rigorous stability proof of a lossless Euler-Bernoulli beam with 
tip mass which is feedback interconnected with a nonlinear spring-damper system and a strictly 
passive nonlinear dynamical system. Such a configuration comes into play if the tip payload is 
interacting with a nonlinear passive environment, if the (nonlinear) dynamics of the torque and 
force actuators are also taken into account, or for a combination of these cases. It is well known 
that the feedback interconnection of passive systems is passive with a storage function that is 
the sum of the storage functions of all subsystems. In the finite-dimensional case, this property 
is advantageously utilized for the controller design where the storage function usually qualifies 
as an appropriate Lyapunov function candidate. For the infinite-dimensional system under 
consideration, the passivity property still ensures that the storage functional is non-increasing 
along classical solutions, however, it is well known that this does not directly entail asymptotic 
stability. In fact, a crucial step in the stability analysis is to prove the precompactness of the 
trajectories. For linear evolution problems this has been reported in many contributions in 
the literature, but when considering nonlinearities this is much more involved. Under rather 
mild conditions on the parameters and functions appearing in the resulting PDF-ODF model 
representing the overall closed-loop system, global-in-time wellposedness is proven by means of 
semigroup theory and the precompactness of the trajectories is shown by deriving uniform-in- 
time bounds on the solution and its time derivatives. With this, asymptotic stability of classical 
solutions can be guaranteed. 


Appendix A. The Operator 


The system (36) is the mild formulation of the evolution problem (jjp)t = Cpyp with yp = 
[it, £ Tip. Thereby Tip := Hq{0,L) x L^(0,L), and 


C 


p ■ 


u 


V 




V 


- U 

. P 


with the domain 

D{Cp) = np:u£ H^iO,L),v£ H^{0,L), 

ku''{L) + kiu'{L) = 0, Au"'(L) - k 2 u{L) = O}. 
The space T-Lp is equipped with the following inner product: 


p L pL 

{yp,yp)p := A / u''u''dx + p / vvdx 
Jo Jo 


(69) 


-I- kiu'{L)u'{L) + k2u{L)u{L). 


The constants ki,k 2 are defined in (1^ and depend, at first glance, on the fixed yo £ ft 


in the proof of Theorem 5.6 Hence, D(Ap) and the above inner product also depend on yg. 


But this does not cause any problems. Anyhow, Step 2 in the proof of Theorem |5.6| shows that 
uo(L) = Uq{L) = 0. Hence, 

We have the following results: 


Kj=K,. 


Lemma A.l. The operator : Tip D{Cp) exists and is a bijection. Furthermore, is 
compact in Tip. 
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Proof. The proof is analogous to the proof of Theorem 3.1 see also Section 4.2 in HU. □ 
Lemma A.2. The operator Cp is skew-adjoint. 

Proof. First we show that Cp is skew-symmetric, i.e. for all y,y& D{Cp) there holds {Cpy., y)p = 

~{y-: ^py)p'- 

pL pL 

{Cpy,y)p = A / v"u"dx — A / u^^vdxKiv'{L)u'{L)K 2 v{L)u{L) 


= A( / vif^ dxv'[L)u'{L) — v{L)u''{L) 


' dx — u'"{L)v{L) u"{L)v'{L^ 


-\- Kiv'{L)u'{L) K 2 v{L)u{L). 


Using the boundary conditions Au"{L) -\- Kiu'{L) = 0 and Au"'{L) — K 2 u{L) = 0 from D{Cp) 
we obtain: 

{Cpy,y)p = A ( vvF^ dx — Kiv'{L)u'{L) — K 2 v[L)u{L) 

Jo 

— A I u"v" da: — K 2 u{L)v{L) — Kiu'{L)v'[L) 

Jo 

d- Kiv'{L)u'{L) -\- K 2 v{L)u{L) 

~ ~ ( 2 /) J^py)p- 

So Cp is skew-symmetric. Furthermore, due to Lemma |A.1| we know that ran£p = Tip. So 
we can apply the Corollary of Theorem VII.3.1 in [3 Uj which proves the skew-adjointness of 
Cp. □ 

Lemma A.3. Cp generates a Co-semigroup of unitary operators in Hp. 

Proof. Since Cp is skew-adjoint, this follows from Stone’s theorem [33l Theorem II.3.24]. □ 
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